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Abstract. In this article we propose a generalisation of the recent work of Gatheral-Jacquier 1111 on 
explicit arbitrage-free parameterisations of implied volatility surfaces. We also discuss extensively the 
notion of arbitrage freeness and Roger Lee's moment formula using the recent analysis by Roper |18| . 
We further exhibit an arbitrage-free volatility surface different from Gatheral's SVI parameterisation. 

1. Introduction 

European option prices are usually quoted in terms of the corresponding implied volatility, and over 
the last decade a large number of papers (both from practitioners and academics) has focused on un- 
derstanding its behaviour and characteristics. The most important directions have been towards (i) 
understanding the behaviour of the implied volatility in a given model (see [2] , [3J [5] , [T^] for instance) 
and (ii) deciphering its behaviour in a model- independent way, as in fTT], [T8] or [JJ5]. These results 
have provided us with a set of tools and methods to check whether a given parameterisation is free of 
arbitrage or not. In particular, given a set of observed data (say European calls and puts for different 
strikes and maturities), it is of fundamental importance to determine a methodology ensuring that both 
interpolation and extrapolation of this data are also arbitrage-free. Such approaches have been carried 
out for instance in [4] , in [7] and in |20j . Several parameterisations of the implied volatility surface have 
now become popular, in particular [13] . [T5] and [9], albeit not ensuring absence of arbitrage. 

In the recent paper Gatheral and Jacquier proposed a new class of SVI implied volatility pa- 
rameterisation, originally proposed in [S]. In particular they provide explicit sufficient and — in a certain 
sense — almost necessary conditions ensuring that such a surface is free of arbitrage. This new class de- 
pends on the maturity and can hence be used to model the whole volatility surface, and not a single 
slice. It also depends on the at-thc- money total implied variance 9t, and on a positive function tp such 
that the total variance w as a function of time-to- maturity t and log-(forward)-moncyness k is given by 
w(k,t) = t SVlp(kip(8 t )), where SVI p is the classical (normalised) SVI parameterisation from fTT], and p 
an asymmetry parameter (essentially playing the role of the correlation between spot and volatility in 
stochastic volatility models). 

In this work, we generalise their framework to volatility surfaces paramcteriscd as w(k, t) = 6t^(kip(9 t )) 
for some (general) function vp. We obtain (Sections [3J and 0]) necessary and sufficient conditions coupling 
the functions "J and (p, so that the corresponding implied volatility surface is free of arbitrage. This 
allows us to obtain (i) the exact set of admissible functions ip in the symmetric (p = 0) SVI case, and (ii) 
a constraint-free parameterisation of Gatheral-Jacquier functions satisfying [IT] Corollary 1]. Along the 
road (Section [5]) we also refine the notion of arbitrage freeness used in fTT] , in connection with the recent 
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analysis by Roper [18]. More precisely, we get a sufficient condition for a real function to represent the 
slice of an arbitrage-free implied volatility surface at a given maturity. We provide a purely analytical 
'ticket proof of Roger Lee moment formula [T7] in our (restricted) case. In the last section (Section [5]), 
we exhibit examples of non-SVI arbitrage-free implied volatility surface after formulating an extension of 
Roper's theorem in implied volatility to convex non-smooth smiles. 

Notations: In this paper, we consider European option prices with maturity t > and strike K > 0, 
written on an underlying stock (St)t>o ■ Without loss of generality we shall always assume that So = 1 and 
that interest rates are null, and hence the log (forward) moneyness reads k := \og(K). We shall denote 
the Black-Scholes value for a European Call option Cbs(^', w), for a strike K and total variance w, and 
more generally C(K,t) for (any) European Call prices with strike K and maturity t. For any k € R 
and t> 0, the corresponding implied volatility is denoted by o~(k,t) and the total variance w is defined 
by w(k,t) := a(k,t) 2 t. With a slight abuse of language (commonly accepted in the finance jargon), we 
shall refer to the two-dimensional map (fc, t) t— > w(k, t) as the (implied) volatility surface. Finally, for two 
functions g and h, we shall say that g(z) ~ h(z) at z = whenever lim z _>o 9( z )/h(z) = !• 

2. Gatheral-Jacquier SVI parameterisation 

Extending the pioneering work of [3], Gatheral and Jacquier introduced in |llj the following parame- 
terisation for the total implied variance w: 



with 9t > for t > and <p is a smooth function from R!j_ to R + and p G (—1,1). The main result in 
their paper (Corollary 5.1) is the following theorem, which provides sufficient conditions for the implied 
volatility surface w to be free of static arbitrage: 

Theorem 2.1. The surface <|2 . 1 1) is 'free of static arbitrage' if the following conditions are satisfied: 



(1) d t 9 t > for allt> 0; 

(2) tp{0) + eLp'{6) > for all 6 > 0; 

(3) tp'(6) < for all 6 > 0; 

(4) Oip(0)(l + \p\) < 4 for all 6 > 0; 

(5) 6ip(9) 2 (l + \p\) < 4 for all 9 > 0. 



A few remarks arc in order here: 

(1) The notion of 'static arbitrage freeness' is defined in [TT| as equivalent to the positivity of the 
numerator and the denominator of Dupire's equation. These properties are necessary, yet the 
right statement is that there exists an inhomogeneous diffusion process associated to the implied 
volatility through Dupire's equation. To this end, the local volatility defined through Dupire's 
equation should have adequate continuity and growth properties, which should be discussed at 
least. An alternative route is to follow Roper conditions in |18j . which add some limiting behaviour 
to Gatheral and Jacquier conditions; we shall come back on this issue in detail in Section [5l 

(2) the conditions in Theorem 12. II are sufficient, but may be far from necessary; 

(3) the full characterisation of the functions tp guaranteeing absence of (static or not) arbitrage in 
the symmetric SVI case p = is left open. 



(2.1) 
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(4) it would be useful to 'parameterise' the set of functions ip satisfying the conditions of Theorcm l2.1[ 
This could lead to easy-to-implement calibration algorithms among the whole admissible class, 
without being tied to a particular family as in [llj . 

In this paper, we try to settle all these points. We also state our results in a slightly more general 
framework, which is not tied to the specific shape of the SVI model, i.e. we consider implied volatility 
surfaces of the form 

(2.2) w(k,t) = 6 t V(kip(9 t )), for all k e R, t > 0. 

We shall make the following assumptions on the parameters: 

Assumption 2.2. 

(i) 9 e C 1 (R* + -> R + ), limt-^o t = and 6 is not constant; 

(ii) ipeC 1 (M* + -> 

(iii) $eC 2 (l^ W + ), *(0) = 1 and * is not constant. 

The time-dependent function 8 models the at-thc-moncy total variance; the assumption on its be- 
haviour at the origin is thus natural. A constant function \t corresponds to deterministic time-dependent 
volatility, a trivial case we rule out here. Likewise, were 9 assumed to be constant, it would be null 
everywhere, which we shall also not consider. The main goal here is to provide sufficient conditions 
on the triplet (9, cp, \&) that will guarantee absence of static arbitrage. Note that when the corre- 
lation \p\ is strictly small than 1, the SVI parameterisation (|2.ip corresponds to the case ^(z) = 
|(1 + pz + y/ z 2 + 2pz + 1). In the sequel, we shall refer to this case as the SVI case. In the follow- 
ing sections [3] and 0] we stick to the definitions of arbitrage freeness as used in [TT| . In Section [5] the 
cleaner formulation of [18] is applied. 

3. Elimination of calendar spread arbitrage 

We first concentrate here on determining (necessary and sufficient) conditions on the triplet (9, <p, ^) 
in order to eliminate calendar spread arbitrage. 

3.1. The first coupling condition. The following definition is now well-known, and further details 
(and a proof) can be found in [TU] or [TT] 

Definition 3.1. A volatility surface w is free of calendar spread arbitrage if and only if dtw(k,t) > 0, 
for all k € R and t> 0. 

The quantity dtw{k, t) is nothing else than the numerator of the local volatility expressed in terms of the 
implied volatility, i.e. Dupire's formula (see QU])- Define now the functions F : R — > R and / : R — > R by 

(3-1) ^) : =4rV> /(«)- 

V(z) ip(u) 

These functions will play a major role in our analysis, and Assumption (2^21 (iii) implies in particular that 
F(z) ~ ^'(0)2; at the origin and F(0) = 0. Note that \1/ and ip can be recovered through the identities 



= exp ^ 



(p{u) = <p{r) exp ' 
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for some arbitrary constant r > 0. The following proposition gives new conditions for absence of calendar 
spread arbitrage. 

Proposition 3.2. (First coupling condition,) Let 0^ := lim^oo 9 t . The surface w in (|2.2[) is free of 
calendar spread arbitrage if and only if the following two conditions hold: 
(i) is not decreasing; 

(ii) \ + F(z)f(u) > for any z eM and u G (O,0oo). 

Proof. By Definition 13.11 the surface denned by (|2.2p is free of calendar spread arbitrage if and only if 

(3.2) d t w(k,t) = 0' t y(z) + 6 t y'(z)ktp'(6 t )9' t > 0, for aU k G R,t > 0, 

where we define z := kip(9 t ). Since $(z) is strictly positive by Assumption 12.21 (iii). the inequality (|3.2p 
is equivalent to Q' t (1 + F(z)f(8 t )) > for all z G R and t > 0, where the functions F and / are defined 
in (|3.ip . For z = we get 0{ > for all t > 0. Otherwise (ii) is necessary and sufficient for the surface 
to be free of calendar spread arbitrage. □ 

Remark 3.3. We do not assume here that 9^ is infinite. In most popular stochastic volatility models with 
or without jumps, 9^ is infinite. In fact, as studied in |16j . as soon as the random family (t^ 1 log(St)j Q 
satisfies a large deviations principle as t tends to infinity with a continuous rate function and speed t , 
then (9oo = oo. Rogers & Tehranchi [19| showed that for a non-negative martingale (St)t>o satisfying 
St > almost surely for all t > 0, the equality 9oo = oo is equivalent to the almost sure equality 
limtyioo St = (where the limit exists by the martingale convergence theorem). However, it may occur 
that < oo. As a corollary of coupling properties of stochastic volatility models, Hobson [T3] provides 
instances where such a phenomenon appears. Simple examples are the SABR |13| model with [3 = 1 and 
models with bubbles. 

Remark 3.4. Condition (ii) in Proposition 13.21 can be stated in a more compact way: 

1 — sup F + sup /_ > and 1 — sup -Fl sup / + > 0, 
where /+ := max(/, 0) and /_ := max(— /, 0). 

Motivated by the celebrated moment formula in [17] , which forces the function \P to be at most linear 
at (plus/minus) infinity, let us propose the following definition: 

Definition 3.5. We say that the function ^ is (strictly) asymptotically linear if and only if ^' has a 
finite (non-zero) limit a+ (resp. a_) when z — > +oo (resp. z — > — oo). 

Using this definition we obtain a necessary condition on the behaviour of the function tp in (|2.2p . 

Proposition 3.6. //vp is strictly asymptotically linear and if there is no calendar spread arbitrage, then 
the map u <— > u<p(u) is non- decreasing. 

Proof. Using Q3.ip . if ^ is strictly asymptotically linear, then \im z ^± 00 z^/'(z)/ 1 ^(z) = \im z ^± 00 F(z) = 
1, so that absence of calendar spread arbitrage implies 1 + f{u) > for any u £ (0,0 oo) by Proposi- 
tion 13.21 (ii) . Since tp is a strictly positive function by Assumption 12.21 (ii) , the proposition follows from 
the definition of the function / in (|3.1|) . □ 
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Note that if lim z _ i .± 00 $>'(z) — then the limit of the function F at (plus or minus) infinity does not 
necessarily exist. Whenever it does, since the function z i-> ^f(z)/z is decreasing as z — > ±00, the limit 
can take any value in (—00, 1). 

3.2. Application to SVI. In the SVI case ([23]) . we have #'(2) = § ^0 + ^ = =^ , so that * is 
strictly asymptotically linear when \p\ < 1 with a + = p + 1 and a_ = p — 1. Therefore Proposition 13.61 
applies, and a necessary condition is that u 1-4 uip(u) is not decreasing. In [11[ Theorem 5.1], this condition 
and the condition that tp is not increasing are shown to be sufficient. In the case of the symmetric SVI 
model (where p = 0), the following corollary relates our conditions and that of [TT] . 



Corollary 3.7. In the symmetric SVI case, the necessary condition of Provosition \376\ is also sufficient 

Proof. Let us consider the the symmetric case p — 0. We compute explicitly 
(3.3) 

1 4- \/l 4- 7 2 1 Z 7 2 

Wz) = +Vi + Z , V(z) = -^^=, *"( z ) = ___i foralUeR. 

Therefore 

F W = . _ and *•(,)_ (~' 2 + 2(l + ^)) = 



vT+7»(l + 1^+72) (l + z 2 ) 3/2 (l + N /TTS 5 ) 

It is then clear that the even function F is strictly increasing on M* + and strictly decreasing on M*_ with 
a global minimum attained at the origin for which F(0) = 0. In light of Rcmark l3.41 we have supF + = 1 
and supF_ = 0. By Proposition ^. 21 there is hence no calendar spread arbitrage if and only if f(u) > —1, 
which is equivalent toiih) u<p(u) being non-decreasing. □ 

4. Elimination of butterfly arbitrage 
We now consider the butterfly arbitrage which, probably not surprisingly, is more subtle to handle. 



4.1. The second coupling condition. We define the operator C acting on C 2,1 functions by 
(4.1) Cw{k,t) := (l - 



kd k w(k,t)\ 2 d k w{k 7 t) 2 (I l\ d 2 kk w{k,t) 



2w(k,t) J 4 \w(k,t) 4, 

for all k £ E and t > 0. The operator does not act on the second component of the function. We however 
keep this notation for clarity. When the function w represents the total implied variance, Cw is nothing 
else that the denominator of the local variance expressed in terms of the implied volatility, i.e. Dupirc's 
formula, see for instance [TUJ Chapter 1]. Recall the following definition from [TTJ Definition 4.1]: 



upi 

B 



Definition 4.1. A volatility surface w is free of butterfly arbitrage if and only if Cw(k,t) > for all 
k e R and t > 0. 

We now reformulate the butterfly arbitrage condition in our setting. For any u > define the set 
(4.2) Z+W := {* 6 > : ^ (^f _ + *t£ > „ 

where we use the convention inf0 = 00. 



^The authors of have informed us that a revised version of the preprint includes an additional condition ensuring 
that the corresponding density integrates to one. Section \E\ below precisely deals with this, so we leave it as it is for now. 
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Proposition 4.2. Second coupling condition, general formulation The surface w given in (|2 .2[) 
is free of butterfly arbitrage if and only if 

2 



1 



s*'(z) 



(uld(u)) 2 < inf ~, —4 , for all u G (0,6L„). 

Au \ Wjz) ^ \ Z ) J ^ 16 

Proof. From we clearly have d k w(k,t) = 6^' (z)tp(6 t ) and dl k w{k,t) = t ^>" {z)tp(6 t ) 2 for all k e 

and t > 0. Therefore 

£w(fc / _ kd kW (k,t)\ 2 d k w(k,tf ( i , i^ , 9f feW (fc,t) 



2w(k,t) J 4 \w(k,t) A J 2 



26> t *(z) / 4 V 6 't*( z ) 4 



2*(z) / 4 V *0) 7 16 

(4.3) = (I-HtttV- li- - W ' (r,2W 



2<5(z) / v ^ v " [46» t V *(*) J 16 

Using the definition of the set Z + (u) in (|4.2p . the proposition follows. □ 

4.2. Asymptotically linear ^P. We now consider the case where ^ is strictly asymptotically linear. 
Define the sets 

(4.4) Z+ := jz G R : (^Jy- - 2¥"(z)) > o| , and Z_ := R \ Z+, 
as well as the constant 

(4.5) Mqo ; = li m uLp(u). 

The following proposition is a reformulation of Proposition 13.61 in the asymptotically linear case. It 
provides sufficient and necessary conditions for the surface (|2.2p to be free of calendar spread arbitrage. 

Proposition 4.3. Assume ^ is strictly asymptotically linear and there is no calendar spread arbitrage. 
There is no butterfly arbitrage if and only if the following two conditions hold: 



ML < inf 



2*(s) 



where is defined in (|4.5[) 
(ii) for any u G (0, 0oo), 

( i _ z*'(z) \ 2 

(utpiu)) 2 < inf — , ^ mz) { 

Moreover Z + is not empty and not bounded from above. 
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Proof. Assume that \& is strictly asymptotically linear and that there is no calendar spread arbitrage. 

Then u i-> uip(u) is non-decreasing on by Proposition l3.6l On Z_ the map u i-» ^ (~¥^) 2\&"(z)J + 

- ^) is clearly also non-decreasing on W + . Therefore (Z + (u)) u> o is a non-decreasing family of sets and 
Z+(0) = Z+. The inequality 

{uip{u)f < 
is equivalent to 



inf 





V 2*0) ) 




J- 

4u 


1 *(z 


2 2*"(z) N 


1 T 16 



for all u e (0, 6>oo) 



Note that in case O 



< inf 
oo this simplifies to 



2*(s) 



2 ez_ n^+^oo) i 



Moo < inf 

zez_ 



- **"(*) 



2 2 



1"(2 



16 



$'(z) *(z) 

Statement (ii) in the proposition is then an immediate consequence of Proposition ^. 2l and the computation 
above. Let us remark also that in this case Z + is not empty. If it is then we can then choose a > such 
that *'(z) > for all z > a. This implies that <3>(z) _1 < 2v&"(z)/*'(z) 2 for all z > a, which in turn 
yields J* *(fe)- 1 d6 < 2 (tf'(a)- 1 - l-'(z)- 1 ). The integral diverges to infinity as z tends to infinity since 
^{z) ~ a + z whereas the right-hand side is bounded by Assumption 12.21 The same argument shows that 
Z-L. is not bounded from above. □ 



Coupling condition (ii) with 



oo, we can also take u to oo, and we obtain the necessary condition 
4 2z 



< inf 



in particular Mqo < 2/ sup{|a + |, 



4.3. Application to symmetric SVI. As in Section l3~2l above, we show that in the symmetric SVI 
case, all our expressions above are easily computed and give rise to simple formulations. Let us define 
the functions Y, A and and A* by 

16uy(y + 1) 



(4.6) 

(4.7) 
(4.8) 



A{y,u) 
Y(u) := 

K ' l-u/4 V \1-uJA) l-u/4 
A*(u) := A(Y(u),u). 



8(y - 2) + uy(y - 1)' 
2 




Of course we only define these functions on their effective domains, the forms of which we omit for clarity. 
The following proposition makes the conditions of Proposition 14.31 explicit in the uncorrclatcd SVI case. 



Proposition 4.4. In the symmetric SVI (|2.1[) case p = 0, there is no butterfly arbitrage if and only if 
uip(u) < A*(w)l {u<4 } + 161 {u >4}, for all u G (0, 6»oo). 
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Proof. For any z G M, using Q3.3P and letting £ x := z 2 and y 2 := \fl + £ 2 , we can compute 

(*'(^)) 2 - 2*(e B )*"«.) = 2 (i/ '~ 2)( ?' + 1)2 ' 

4 

1 z *'^ = l( 1 + ± 



2*(&) 2 V 2/. 

where by ^'(^z) we really mean 9 u v E'(u)| M= ^ . From the first equation, the set Z + defined in (|4.4p is 
exactly {z G K : y z > 2} = R\ [— V3, V3]- For any fixed u, the function appearing on the right-hand side 
of Proposition 14.31 (ii) can be computed explicitly and is equal (after simplifications) to A(y Zl u) given 
in (|4.3p . In particular A(2, u) = 48 and lim^oo A(y, u) = 16. For any u > 0, we have 

OyA(y,u) = 1, 

(8y — 16 + y 2 u — yit) 

where := (l - f ) y 2 - 4y - 2. When u>A,B u is concave on (2, oo) with JB U (2) = -(6 + u) < 0, 

and hence the map y t— > A(y, u) is decreasing on (2, oo) and its infimum is equal to lim^oo A(y, u) = 16. 
For u G [0,4), the strict convexity of B u and the inequality B u (2) = —(6 + u) < implies that the 
equation B u (y) = has a unique solution in (2, oo), which in fact is equal to Y(u) given in (|4.7[) . Then 
the map y A(y,u) is decreasing on (2,Y(u)) and increasing on (y(u),oo). Its infimum is attained at 
Y(u) and is equal to A*(u) defined in (f478|) . □ 

Remark 4.5. In jTT] , the authors prove that the two conditions (altogether) u(p(u) < 4 and uip(u) 2 < 4 
(for all u > 0) arc sufficient to prevent butterfly arbitrage in the uncorrelatcd p — case. These two 
conditions can be combined to obtain (uip(u)) 2 < 16 min(l, (^(m) _1 ). A tedious yet straightforward 
computation shows that A* is increasing on [0,4) and maps this interval to [0, 16). Notwithstanding the 
fact that our condition is necessary and sufficient, it is then clear that 

(i) for u > 4, it is also stronger than the one in |llj whenever tp(u) < 1; 

(ii) for u < 4 (which corresponds to most practically relevant cases) it is also stronger whenever 
16/ <p(u) < A*(u). 

In particular, item (ii) could be used as a (moving) sufficient and necessary lower bound condition for 
the function (p on [0,4). 

5. Generalised Gatheral-Jacquier surfaces and absence of arbitrage 

5.1. Roper Large-Moneyness Behaviour condition. As recalled in Definitions 13 . 1 1 and 14.11 Gathcral 
and Jacquier |llj define the absence of calendar spread arbitrage and of butterfly arbitrage by the non- 
negativity of the numerator and denominator of Dupire's equation in implied volatility. These properties 
are necessary, yet not sufficient, as shown by Roper [18]. More precisely, the additional condition from p~8| 
Theorem 2.9], namely the Large-Moneyness Behaviour (LMB), reads 

(5.1) lim w(k, t)) = — oo, 

fc/"+oo 

with the usual notation from the Black-Scholes formula d±(k,w) := —k/y/w ± s/w/2. Note that this 
condition is equivalent to call option prices tending to zero as the strike tends to (positive) infinity, as 
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proved in [19l Theorem 5.3]. As mentioned in the footnote on page[5l this condition has just been added to 
Definition ^. ll in a revised version of [TT|. The following lemma however shows that some other asymptotic 
behaviours of d + and d- hold in great generality. This was proved by Rogers and Tchranchi [T§] in a 
general framework, and we include here a short self-contained proof. 

Lemma 5.1. Let w be any positive real function. Then 
(i) lim fe /. +00 d-(k, w(k)) = -oo; 
(ii) lirrifc^-oo d + (k, w(k)) = +oo. 

Proof. The arithmetic mean-gcomctric mean inequality implies 



-*_(*,«;(*)) = -jt= + > V2k, forfc>0, 

V w ( k ) _z 

d+{k, w{k)) = + ^ > V^2k, for k < 0. 

V w ( k ) 2 

In the first line, the term on the right-hand side tends to +00 as k tends to +00, proving (i). Statement (ii) 
follows by the same argument. □ 

The following proposition bridges the gap between our framework and [18\ Theorem 2.9]: 

Proposition 5.2. Assume in the setting of Proposition that the inequality in Condition (ii) is strict. 
Then the corresponding implied volatility surface is free of static arbitrage. 

Proof. In our setting (with being strictly asymptotically linear) \mik^+ cx> (w(k,t)/k) = 6 t tp{0t)ct+. so 
that all we need to prove is that 6 t <f(0 t ) < 2/a+. For any z S Z + (defined in (|4.4|Q . note that 

2 / T ,, , \ 2 



1 _ *fr'(z) 
1 2*(z) 



< 



f 1 _ rf'W V 

V 2*( 2 ) ) 



j_ ( *>(z) 2 _ 2^/(^0 

Choosing a sequence in Z + diverging to infinity, we have (9t^p{6t)) 2 < 4/a^_ and the proposition follows. 

□ 

5.2. Small-moneyness behaviour, no mass at zero, and symmetry. The missing statements in 
Lemma I5TT1 arc the LMB Condition (|5.1[) and the Small-Moneyness Behaviour (SMB): 

lim d-(k,v(k)) = +00. 

To investigate further, let us remark that the framework developed in |18j encompasses situations where 
the underlying stock price can be null with positive probability. This can indeed be useful to model the 
probability of default of the under Iyer. Computations similar in spirit to |18j show that the marginal law 
of the stock price at some fixed time t > has no mass at zero if and only if 

lim d K C(K,t) = -1. 

which is a statement about a 'small-moneyness' behaviour. This can be fully recast in terms of implied 
volatility, and the above missing conditions then come naturally into play: 



Proposition 5.3. (Symmetry under small-moneyness behaviour) 
Let a C 2 real function v satisfy 
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(I) v(k) > and Cv{k) > for all k e M; 
(II) linifeN^-oo d— (k, v(k)) = +oo (SMB Condition); 
(III) limfcy! +00 d.\-(k, v(k)) = — oo (LMB Condition). 
Define two functions p- and p + by k i— > p±(k) := (2i:v(k)) cxp (— \d\ (fc, Cv(k). Then 

(1) p+ and p- define two densities of probability measures on K with respect to the Lebesgue measure, 
i.e. j R p^(k)dk = J R p+(k)dk = 1. 

(2) p+{k) = c fe p-(fc), so that X!^ c k p_(k)dk = c~ k p+{k)dk = 1. 

(3) p- is the density of probability associated to Call option prices with implied volatility v, in the 
sense that p-(k) = dj (K CBS (K, v(log(K)))\ K k , and k i— > p + (—k) is the density of probability 
associated to Call option prices with implied volatility k t— > w(k) := v(—k). 

Remark 5.4. Note that the strict positivity of the function v in Assumption (I) ensures that the support 
of the underlying distribution is the whole real line. One could bypass this assumption by considering 
finite support as in [TS] . In the latter case — slightly more general case — the statements and proofs would 
be very analogous but much more notationally inconvenient to write down. 

Proof. The functions p_ and p+ are clearly well-defined and non-negative. Consider first p_. It is readily 
seen that the function D{k) = Sr-Cbs {K, v(\og(K)))\ K=ck is a primitive of p_. We now proceed to prove 
that p_ is indeed a density. Let N denote the cumulative distribution function of the standard Gaussian 
distribution. An explicit computation yields (the reverse one can be found in [111 Lemma 4.1]) 

d K C BS {K,v(log(K))) = d K d + -^= G - d +' 2 -N{d-) - Kd K d--^=e- d2 -' 2 

V27T V27T 



l — e ~ d \n (d K d+ - Kd K d-c d +^- v l' 2 \ -M{d 



V2tt 

= -^e" d + /2 {d K d+ - died-) — J\f(d-), 

V 2-7T 

where all the functions are evaluated at (K, log (K)), and where we have used the equality cL = d+ — t/v. 
Evaluating the right-hand side at K = e k , and using the fact that — k — \d\ = — \d?_, we obtain 

D(k) = f<*> ex P (y-^ 2 ) 



2y/2nv(k) 



Therefore if 



(5.2) lixn ex P f- rf -^ fc )) 2 )=0, 

V ' k^±^2^vjk) \ 2 



then 



p-(k)dk= lim Af{d_{k,v{k))) - lira Af(d-(k,v(k))) = 1, 



where we have used Assumption (II): linifc^-oo d— (k, v(k)) = +oo. Let us first deal with the case where k 
tends to (positive) infinity. From Lemma [5.11 cxp (— \d- (k, f(fc)) 2 ) tends to zero. The key point is to 
observe that D is the primitive of a non-negative function, therefore it is non-decreasing and has a 
(generalised) limit L e [—00,00] as k /• +00. Since M (d- (k, v(k))) converges to zero by Lemma [54] (ii) , 
we deduce that ^-^Ll£L= cxp (— |gL (k, v(k)) 2 ) tends to L. From ]T9j Proof of Theorem 5.3], the inequality 

v' < \j2v jk holds for any k > so that L is necessarily non-positive. Assume that L is negative (possibly 
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infinite) . Since v ^ = dk \j v(k), if L is negative then ^/v is eventually decreasing. Since it is bounded 




from below by 0, there is a sequence (fc n )n>o going to infinity such that dky/ v(k n ) converges to zero by 



the mean value theorem. Therefore L must be 0. 

Let us now consider the case where k tends to negative infinity. Using similar arguments, the quantity 



k small enough. Since v is positive, this implies that v(k) > e for some e > and k small enough. In 
particular l/v(k) is bounded. Since v'(k) > —4 for all k £ R by [T9J Theorem 5.1], then for k small 
enough, the inequalities —4 < v'(k) < hold, and the term outside the exponential in (|5.2p is bounded. 
Since the exponential converges to zero by Lemma 15.11 (ii), we obtain M = 0. If M > 0, then v is 
increasing for k small enough. We conclude as above by the mean value theorem since y/v is increasing 
and bounded from below. Therefore M = and the limit ()5.2[) holds. 

So far we have proved that p_ is the density of probability associated to Call option prices with implied 
volatility k v(k). Consider now the function w(k) = (—k). Then for all k £ R, dkiv(k) = —dkv(-k), 
and it follows by inspection that Cw{k) = Cv{—k) > 0. Consider the function p_ associated to w, i.e. 



Now d-(k,w(k)) = fc) _ 2 = -C M — fc> u (~ 80 tnat P-W = P+(—k). In order for 

p_ to be a genuine density, we need to check conditions symmetric to those ensuring that p_ is a 
density. The condition symmetric to the SMB assumption (II) is precisely Condition (i) in Lemma |5.1[ 
and the condition symmetric to the Lemma 15.11 (ii) is precisely the LMB assumption (III) . Therefore 
k i— > p+(— k) is also a density, associated to Call option prices with implied volatility w. Finally the 
identity p+(k) = c k p-(k) follows immediately from the equality — k — hd?L = — ^d 2 _. □ 

The previous proposition has been intentionally stated in a maturity-free way: it is indeed a purely 
'marginal' statement. A natural question is then to determine whether such a function v, satisfying the 
assumptions of Proposition 15.31 can represent the total implied variance smile at time 1 associated to 
some martingale (issued from 1 at time 0)? The answer is indeed positive and this can be proved as 
follows. Consider the natural filtration B of a standard (one-dimensional) Brownian motion, (B t )t>o- 
Let P be the cumulated distribution function associated to p- characterised in Proposition 15. 3[ and let 
Af be the Gaussian cumulative distribution function. Then the random variable X := P~ 1 (Af(Bi)) has 
law P, and E(X) = 1. Set now M s := E(X\B S ), then M is a martingale issued from 1. Note that M 
is even a Brownian martingale and therefore a continuous martingale. The associated Call option prices 
E[(M S — K) + ] uniquely determine a total implied variance surface (t, k) <— > w(k, t) such that v = w(l, ■). 

5.3. A simplified version of Lee moment formula in the strictly asymptotically linear case. 

So far we have stressed that, following Roper or the variant in Proposition l5.3l (no mass at the origin), the 
positivity of the operator £ in (|4.1[) guarantees the existence of a martingale explaining market prices. 
As a consequence, the celebrated moment formula [17] holds. We show in this subsection that, at least 
in the strictly asymptotically linear case, it can be derived in a purely analytic fashion. 

Proposition 5.5. Consider a C 2 function v satisfying the following conditions: 




tends to M £ [—00,00]. Assume that M < 0. Then v is decreasing for 




for all k £ R. 



(1) v(k) > and Cv{k) > for all k £ R; 



12 



GAOYUE GUO, ANTOINE JACQUIER, CLAUDE MARTINI, AND LEO NEUFCOURT 



(2) \\m k s +00 v'{k)=a£ [0,2); 

(3) lim fe ^ +ooW "(fc) =0. 

Let X be a random variable with density p, where p is the probability density associated to v by Proposi- 
tion\5^ Then X satisfies E (X) = 1 and sup {to > : E (X 1+m ) < oo} = - - 1 + a) . 

Proof. Condition (1) clearly implies Condition (I) in Proposition 15.31 and Conditions (2) and (3) imply 
the SMB and LMB limits (II) and (III) in Proposition 15.31 Therefore by Proposition 15 . 31 the centred 
probability density p is well defined on R. For m e I, we therefore have e^ 1+m ^ fc p(fc) = f{k)e~ 9 ^ k \ where 

f(k) = (2irv(k))- 1/2 £v(k) and g(k) = \ (J^- + ^ + - (1 + m)k. 

As k tends to infinity, straightforward computations show that 

f{k)~ 4 7" 2 -d lim ^ = (a~2) 2 -8ma 

16V27rafc fc/-+oo k 8a a 

Since P m is a second-order strictly convex polynomial with P m (0) > 0, the function k i-> /(/c)e^ 1+m ^ fc p(fc) 

is intcgrable as long as P m (a) > 0, i.e. a < 2 — 4 (v^ 



'm + to — to j , or 
a 1 1 

TO < 1 . 

8 2 2a 



In other words, we have proved that: 



sup {m > : E (X 1+m ) < oo} = i - 1 



Q 



□ 



5.4. A variant of Roper's theorem for non-smooth implied volatilities. Let us observe that the 
price formulation of arbitrage freeness in [TSJ Theorem 2.1] is minimal in the sense that the regularity 
conditions on the Call option prices are necessary and sufficient: to be convex in the strike direction and 
non-decreasing in the maturity direction. The implied volatility formulation |18| Theorem 2.9, condition 
IV. 1] however, assumes that the total variance is twice differcntiablc in the strike direction. This regularity 
is certainly not required; in fact, the author jTS] Theorem 2.9] proves the latter by checking the necessary 
assumptions of Theorem 2.1 on the properties of the Call option prices defined by the formula 

(5.3) C(c k ,t) = C B s(c fe ,w(M)), for all k G R,t > 0, 

where Cbs represents the Black-Scholes call price. More precisely, Roper uses the regularity assumption 
in k of w in order to define pointwise the second derivative of C with respect to the strike. He then 
proves that the latter is positive, henceforth obtaining the convexity of the price with respect to the 
strike (Assumption A.l of Theorem 2.1). It turns out that the same result can be obtained without 
this regularity assumption: assume that for any t, the function k \-> w(k, t) is continuous and almost 
everywhere differentiable. Since dkW is defined almost everywhere, and since the term in dkkW is linear, Cw 
can be defined as a distribution and we replace the assumption Cw(k, t) > everywhere by Cw(k, t) > 
in the distribution sense. In order to do so on R* — saving us from dealing with the boundary behaviour 
at the origin — we can assume additionally the SMB Condition (Proposition 15.31 (II)) and work on R^_. 
To conclude, we need to prove that the Call options C defined in (|5.3[) are convex function in K = c k . 
Since a continuous function is convex if and only if its second derivative as a distribution is a a positive 
distribution, the only remaining point to check is that the second derivative as a distribution of C is Cw. 
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But this the same computation already carried out in |18| . and the result hence follows. Let us finally 
note that our assumptions on w are indeed minimal: conversely, if we start from option prices convex 
in K, their first derivative are defined almost everywhere, and so are those of w (in K or k) since the 
Black-Scholcs mapping in total variance is smooth. 

6. The quest of ^ alternative to SVI 

In order to find examples of pairs ip, with "J different from the SVI parameterisation (|2.ip . observe 
first that the second coupling condition (Proposition 14.2] ) is more geared towards finding out ip given \E' 
than the other way round. We first start with a partial result in the other direction. 

6.1. An effective bound for "J. Wc consider the quasi-linear case lim w'(k) = a+ > 0, when 6^ = oo. 

feyoo 



Proposition 6.1. Let w be the generalised SVI surface (|2 .2[) . // the surface is free of static arbitrage, 
then there exists z+ > and k > such that for all z > z + the upper bound 



V ' - M„ V Moo 



holds, where is defined in 

Proof. In the generalised SVI case (|2.2[) . the function "J is strictly asymptotically linear (see Definition l3.5p 
with liuiz^ao fy'(z) = a + > 0, and 6^ = oo. From subsection 14.21 the condition < 

holds for all z € Z + (9 ao ) = R. Since lim Z /»oo f^T^y — l^fyl = ;^r> we can define 

z+ :=i„f{ 2eR+ :i § f( 5 i ;J -^ J )>o}<oo, 



4 2z 
WJz) ~ W(zj 



and therefore 



(6.1) _t__J^_> M _ forallz>z + . 

Note in passing that Moo < 2/a+. Define u + := fy(z + ). Since the (continuous) function 'J is increasing on 
[z + , +oo), we can define its inverse g : [u + , +oo) — > [z+, +oo). Using the equality cxp = 
we can rewrite (|6.1[) as 

4 2z fl (u) Moo 

> Alee .g (u) - — — > — — 



~*H-(-jC£))**W? 

, x / ti -J- , . JV/qq / / \ 

9(u)\ — -g{u + ) > —^y/u^(*/u- y/u^) 

y u 2 

> g(u+)\ + —^-Vu(Vu- y/u+) . 

V u+ 2 

where all the inequalities on the right-hand side are considered for it > w + . The third line- is obtained by 
integration between w + and u on both sides of second line. Let K\ := Moo/2 and K s := u + g{u + ) — 
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iMooy^uZjl. We then obtain the condition 

(6.2) g(u) > K s y/u~ + K t u, 



for all u > it 



+ ■ 



Note that K s remains non-negative if we increase z + or decrease M^; indeed lim z yi 00 (2z/^E'(2)) = 2/a+, 
so that the condition M m < 2/a+ is equivalent to < 2z/^>(z) — 2g(u)/u. Finally let us translate 
condition (|6.2p into conditions on \E'. Fix u > u + and denote z := <?(u), then (z — Kiu) 2 > K 2 u, which is 
equivalent to Kfu 2 — (K 2 + 2Kiz)u + z 2 > 0. The discriminant is equal to K 2 (K 2 + 4zKi) and is clearly 
non-negative. Condition (|6.2[) is therefore equivalent to 

X 2 + 2K lZ - K s ^K 2 + 2K lZ K 2 + 2K lZ + K s ^/K 2 + 2K l z 



2J? 



Given z — A;u > (equivalcntly ^(z) < z/Ki) we obtain 'S(z) < k 2 + Az — k\/k 2 + Az, where k 
A' s / (>/2ifi) and A := Kf 1 . 



□ 



6.2. An example of non-SVI function. We now provide a triplet (8,(p,^>), different from the SVI 
form (|2.1[) . which characterises an arbitrage-free volatility surface via (|2.2[) . Define the function # by 
9t = t and the functions ip and \& by 



<^>(u) :={ u ' lf " > ° and *(«) := \z\ + y/l + \z\) , for all z e R. 

( 1, if u = 2 

A few remarks arc in order: 

• the function tp is continuous on R + ; 

• ^ = oo; 

• the map u i— > utp(u) is increasing and its limit is = 1; 

• the function VP — directly inspired from the computations in Section I6.H — is symmetric and con- 
tinuous on R. It is also C°° on K*, and strictly asymptotically linear. Its derivative has a 
discontinuity at the origin, but the results in Section 15.11 allow it. 

With these functions, the total implied variance (|2.2[) reads 

Vt 



w(k,t) = k(i-e- t ) + ^(yi+^/^r^y+tj , 

Proposition 6.2. The surface w is free of static arbitrage. 



for all k 6 



Proof. The function / defined in (|3.1[) therefore reads f(u)= — and is decreasing from to — 1. 

Regarding the function F, it is clearly continuous, increasing from to 1 and 

IzlfvTTR+i) 

F(z) = \ , ' forallzeR*, 

2 v /ITH(2|z| + i + VTTR) 

with F(0) = 0. By Proposition 13.21 straightforward computations then show that the volatility surface 
is free of calendar-spread arbitrage. Now, for any z > 0, we have 

V{z) 2 _ = (16z + 19)yTT^ + 12z + 10 

*(z) (Z> 16(l + z) 3 / 2 *(z) 



GENERALISED ARBITRAGE-FREE SVI VOLATILITY SURFACES 



15 



which is a decreasing function of z with limit equal to zero. Therefore Z + = M, and for any u £ R!j_, 
Z + (u) = R. Let us check that the generalised SVI surface w parameterised by the previous triplet (9, <p, ^) 
satisfies Cw > as a distribution. Indeed we only checked that {Cw} > as a function defined everywhere 
except at the origin (where as usual in distribution notations {Cw} is a function defined where w" is 
defined). Here = {'&"} + §#o (where So stands for the Dirac mass at the origin), so that Cw = 
{Cw} + 5(9(f>(9)) 2 So, which is positive since {Cw} > 0. Finally, 

4 2z 4(z + l) 3 / 2 + 3z + 4 

W(zj ~ ~ (4 v / ^Tl + 1) (2z + 1 + 7z + f) 

decreases to 2 > M^. Since ^'(z) 2 > 1, the condition (uip(u)) 2 < 4 is sufficient for absence of butterfly 
arbitrage. □ 




Figure 1. Plot of the differential operator C (left) defined in (|4.1|) in the non-SVI case 
of Section [621 an d of the density at time t = 1 (right). Note that Cw(k,t) remains 
positive for all k £ K and t £ K+ as expected. 



6.3. An example of non-SVI parametric family. We propose a new triplet (0, <p, ^) characterising 
an arbitrage-free volatility surface via (|2.2[) . Define the function 9 by 8t = t and 

f l-e~" 

¥>(«):= ^ " u ' ltU>U and 9 v {z) := (1 + \z\ v ) 1/v , for zel, 

la, if u = 

where ^ is a real number in (l,oo) and a € (0,a) with a ps 1.33. Note that when = 2, modulo a 
constant, the function ^2 corresponds to SVI. 



Remark 6.3. We could in principle let a depend on v. The reason for the construction above is that 
we want to show that the corresponding implied volatility surface is free of static arbitrage for all v > 1. 
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The same remarks as in the example in Section [6.21 hold: ip is continuous on R + , 9^ = oo, u i— > utp(u) 
is increasing to Moo = a and \I/„ is symmetric and continuous on R. It is also C°° on R*, C 1 on R, and 
strictly asymptotically linear. With these functions, the total implied variance (|2.2|) reads 



w(k, t) 



1 



(1 



a v \k\ v 



Proposition 6.4. TTie surface w is free of static arbitrage. 



for all k G R, t > 0. 



Proof. The function / defined in (|3.1j) therefore reads /(it) = — -, with /(0) =0 and is strictly 
decreasing from to —1. Regarding the function F, it is clearly continuous, increasing from to 1 and 
F(z) = \z\ v I (1 + \z\ v ) for all z G R*, with F(0) = 0. Since 9. is increasing and 1 + f(u)F(z) > 



for all (u,z) € R^j_ x R, the first coupling conditions in Proposition 13.21 are satisfied, and the volatility 
surface is free of calendar spread arbitrage. We now need to check the second coupling condition, namely 
Proposition 14.21 For v > 2, "J is C 2 , and we can indeed apply Proposition ^. 21 Since \I/ is asymptotically 
linear, we can alternatively check Proposition [473] The equality 

_ KM 2 



(6.3) 



2vi(z) = (i+\zn 1/u - 2 \zr 2 < 



\z\ v -2{v-l)) 



*„(*) 

holds for all z G R* and hence the sets Z + and Z_ defined in (|4.4p are equal to Z_ = [zl , 
Z + = R \ [zl, z+], where z* ± := ±[2(i/ — l)] 1 ^- The two conditions in Proposition 14.31 read 

2 



and 



inf 



1 



z-a' v {z) 

2*„(*) 
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(6.4) 



(u<p(u)) < inf 



1 - 



2*„(z) 



for any u G (0, 00). 



4« »W 1 16 

From the proof of Proposition ^. 3[ we know that when 9 00 = o°j the first condition simplifies to 

4 2z 



(6.5) 



Now, immediate computations yield 



< inf 



2- 



2z 2 + 



\A V (\ + \z\») 1,v 



which, as a function of z is defined on R*, is strictly increasing on W_ and strictly decreasing on Rl. 
Therefore, is infimum z„ over the interval [zl, z* + ] is precisely attained at z± (by symmetry) and is equal 
to 4^(2^ — 2)( x ~ v '/ v (2v — \)~ x l v . Since by construction = a, Inequality (|6.5[) is therefore equivalent 
to a < z„. This inequality is clearly not true for any v > 1 and a > 0; however straightforward 
considerations show that there exists a unique f* > 1 such that the map (l,oo) 3 v t— > z v is strictly 
increasing on (1, v*) and strictly decreasing on (y* , 00) with zi = 4 and liia z ^ +aa z v = 2. Therefore the 
inequality a < z v is satisfied for all v > 1 if and only if a < 2. 

We now check the second inequality (|6.4j) above. Straightforward computations show that 



1 - 



2tf„(z) 



a+i*r 

1 + 
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which increases on R_ from 1/4 to 1 and decreases on R + from 1 to 1/4. The map z i-> ^' u {z) 2 is 
decreasing on R_, increasing on R + and maps the real line to the interval (0,1). Therefore, for any 
u > 0, and any z € Z+, we have 

y 1 - gggyj i i = i 

with defined in (|6.3|) . Now a quick look at the function shows that it is bounded above by 
<& v {zt) G (0,1), with z* := [v(v - 1) - 2 + y/v{v - 1)(^ 2 + Zi> - 2)] 1 /" '. Define the function 5q by 
g a {u) = (ity(w)) 2 (- + ij. There exists a unique u* « 1.87 such that is strictly increasing on (0,u*) 
and strictly decreasing on (u* , oo) with g a (u*) = gi(w*)a 2 . Setting a := gi(u*)~ 1 / 2 w 1.33, the inequality 
ga(u) < 1 is clearly satisfied for any u > and all a £ (0,a). To conclude, note that for 1 < i> < 2, 
the second derivative has a mass at the origin, but ^ „ is convex which implies that this mass is positive 
and that £w(k,t) > in the distributional sense following Section 11x41 Therefore the implied volatility 
surface is free of static arbitrage and the proposition follows. □ 





FIGURE 2. Plot of the differential operator C (left) defined in (|4. 1 [) in the non-SVI case 
of Section 16731 and of the density at time t = 1 (right), with v = 3.5 and a = 1. Note 
that Cw(k,t) remains positive for all k £ R and t £ R+ as expected. Note further that 
in this example, the density does not have a spike at the origin. 



7. Conclusion 

Adopting a slightly more abstract point of view than [TTj , we formulate necessary and sufficient condi- 
tions for a wide family of implied volatility surfaces to be free of arbitrage. This allows us to settle some 
questions raised in [TT|. In particular we characterise the exact set of functions ip for which butterfly 
arbitrage cannot occur in the symmetric SVI model. We precise the notion of arbitrage freeness use in |11) 
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in relation with Roper conditions in |18j , and complement those conditions in order to deal with marginal 
probabilities not weighing zero and non-smooth volatility functions. As an illustration, we provide two 
examples of smooth non-SVI arbitrage-free volatility surfaces. 
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